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PREFACE 


I'  the  devoloiHncnt  of  this  signal  structure  specification  grew  out 

^  of  discussions  between  the  author  and  Major  Nick  Aloxandrow,  then  of  the 

Air  Force  Nuclear  Criteria  Croup  Secretariate.  I  would  like  to  acknowledge 
I  his  contributions  and  the  assistance  of  Dr.  Dennis  Knopp,  Mission  Research 

I  Cor|)oration  and  Dr.  Clifford  I’rettie,  Berkeley  Research  Associates  in 

r  preparation  of  this  report. 
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A  TRANS- IONOSPHERIC,SIGNAL  SPECIFICATION 
FOR  SATELLITE  C^  APPLICATIONS 


1 .  INTROQUCTION 

Proiier  design  of  radio  frequency  systews  that  must  operate 
through  structured  or  varying  plasmas  require  an  accurate  and  easily 
imi>le»entahte  signal  specification.  The  specification  should  be  a 
"reasonable  worst  case"  in  that  any  system  which  operates  with  the 
specified  signal  can  also  operate  ujuler  any  likely  disturbed  condition. 

The  following  specification  is  the  result  of  a  dedicated 
program  s^jonsoved  by  the  Defense  Nuclear  Agency  atjd  the  Air  Force 
Weapons  Laboratory.  Other  ma jar  contributors  were  the  Naval  Research 
Laboratory  and  the  Air  Force  Ceophysics  Laboratory,  the  program  ineU*ded 
investigations  in  the  mechanisms  that  cause  structured  plasmas,  in  the 
propagation  of  electromagnetic  signals  through  structured  plasmas,  and 
in  the  effects  of  scintillated  or  otherwise  distorted  signals  on  typical 
satellite  C  links.  Field  experiments  and  measurement  programs  were 
executed  to  verify  theoretical  plasma  structure  and  radio  propagation 
predictions  ana  to  characterise  the  morphology  of  the  natural  ionosphere. 
Illicit  in  the  resulting  specification  are  several  considerations  and 
asstpptions.  Rather  than  discuss  all  of  the  supporting  logic  leading  to 
this  specification,  a  selected  bibliography  is  included.  Qther  st^iporting 
informatioft  is  available  tn  the  classified  literature.  It  is  noted, 
howev-er.  that  this  specification  should  be  adequate  for  all  satellite 
and  data  link  applications. 


? 

{ 
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2. 


SIGNAL  EFFECTS  FOANALISM 


llie  effects  of  o  disturbed  ionospheric  channel  are  i^epr'^sentod 
by  the  channel  impulse  response  function. 


U(t) 


CO 

J dt  h(t,T)  S(t-T) 


CD 


S(t)  3  transmitted  signal 

R(t)  =  received  signal  (comj^lex  number  representation) 
h(t»r)  a  channel  impulse  response  function 


Tive  specification  chatmel  impulse  response  function  is 


CO 


where  f^  «  carrier  fre^tuency  £hs) 

e  »  light  speed  CSxl©  a/sec) 
sCt)  »  propagation  path  length  Ca) 
r.  »  classical  electron  radius  C2.S-slO  a) 

9 

Sift)  »  total  electron  content  ia”*) 


Far  aiher  than  circular  potarisatioa.  Equation  I  should  he  applied  to  each 
polarisation  state  attd  NCt)  used  to  calculate  the  Faraday  rotation  effects. 
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(3) 


A(t)  »  cxp( 


•0.1iS!C^)i 


oxp(-e^/2o^) 


«  absorption  (dU) 

G(9)  »  antenna  gain  function 

2  2 
cf.  *  energy  angle  of  arrival  variance  (rad  ) 


HJt,4f} 


M  tim 

H 


jdf  h  (f,T)oxp(i2»AfT*i2itft) 


(4) 


Equation  5  assumes  t(utt  the  transmitter  is  on  the  antenna  axis  and  that 
the  antetum  pattern  is  cylindrlcally  symmetric  about  that  axis.  If  the 
threat  has  amplitude  fluctuations  (always  assumed  Rayleigh),  then  h  Cf,T) 
is  a  tero  mean  nortsally  distributed  random  variable  with  an  autocovariance 
defined  by 

h*Cf.^)h^'Cf''rtO  »  *(f-f*)(5(T-ir‘)r2(f.t)  (Sa) 

»  0  (Jb) 

where  is  the  generalised  power  spectrum  and  (!(x-x*)  is  the  Pirae 

delta  function  defined  for  any  ftmetion  f(x)  by 

dx  «(X-X*)fCx)  9  fCx*) 

The  generalised  power  spectrum  is  parameter' ;ed  by  the  scintillated 
signal  decorrelation  time,  and  f  .  the  frequency  seleetiv‘©  batwhsidth. 

Ihere  are  two  forms  for  the  generalised  power  spectrt®  depending  on  the 
product  f^T  where  T  is  the  ainiaua  s>mbol  period.  For  fgT^l 


r2(f.t) 


a 


(6) 


1.864T  5(t) 

0 

2  2T7 
l48.S72T^r  “ 

0  J 


This  represents  the  flat  fade  condition  with  respect  to  the  s/mbol  period, 

T.  for  f  Til 
a 


«  1 

• 

Idx  exp 

4,2 
-X  -2x 

J 

mm 

4  '  Cj  ' 

L 

m 

where  f  »  f^(l*Cj)** 

s  delay  parameter  (=0.25) 

Par  both  equation  6  and  7 


(&) 


iquatien  7  provides  for  seint illation  random  delay  ar,  ©qnivalently. 
frequency  selective  effects. 

If  the  seifttillatioa  threat  consist*  of  phase  effects  only  or 
if  there  are  no  scintillations  at  all*  thee 


S 


(9) 


-  5(05(t) 

If  there  arc  ncintiltations  of  any  kind,  then  there  a  randoa 
time  varying  comjKjnotit  in  N(t)  in  addition  to  the  very  large  s  r  ie 
(slow)  variations. 

The  total  electron  content  is 


N(t) 


df  g(f)exp(-i2itft) 


(Id) 


where  Nj(t)  is  the  large  scale  (slow)  co^jottent  and  g(n  is  a  tero  mean 
normally  distributed  variable  representing  the  random  component.  The 
autocovarlance  of  gif)  is 


whei'e 


f<f 

»  *  *  *t: 


B(f}  ?  g*('f) 


3  Ssyleigh  fretiueacy 


(ti) 


k'ith  the  spoelfiestiea  of  it  M(c}.  ?  .  f„.  f  .  ami  ;(c). 

A  ^  ^  K  6 

we  abtaifi  a  complete  deseriptien  of  the  disturbed  si^t  as  deiersiaed  by 
the  prepagatisa  eaviroRsent  and  systs-m  geometry.  £(&)  is  supplied  by  the 
user  iadepeadeat  of  the  eaviraaeeat  sad  geometry. 
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3. 


SIGNAL  SPECIFICATION  DEFINITION 


In  principle,  a  complete  signal  structure  specification  would 
include  all  of  the  time/ space  variations  possible  for  the  parameters  and 
functions  in  Equation  2.  In  practice,  it  is  sometimes  sufficient  to 
specify  extremum  values  or  simple  functions  for  the  propagation  quantities 
or  their  derivatives  accompanied  by  application  rules.  The  following  table 
lists  a  minimum  set  of  specification  parameters. 

TABLE  1 .  SIGNAL  SPECIFICATION  PARAMETERS 


I 

t 

i 

{ 

i 

h 


I 


SPECIFICATION  PARAMETERS 


Maximum  Values 
Absorption, 

2 

Energy  Angle  of  Arrival  Variance, 
Transraitter/Receiver  Vehicle 

d” 

Dynamics,  — -  ,  n^OjS 
dt" 

Signal  Decorrelation  Time, 


Rayleigh  Frequency,  f^ 

K 

Minimum  Values 

Signal  Decorrelation  Time, 
Frequency  Selective  Bandwidth,  f^ 
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Two  nearly  univot'sal  application  rules  apply  to  and  For 

where  both  a  minimum  and  maximum  is  specified,  the  system  must  operate 

over  all  inteimediatc  values.  Similarly,  when  a  minimum  is 

specified,  the  system  must  handle  all  f^  from  the  minimum  to  the  carrier 

frequency.  Exorcising  those  ranges  are  necessary  because  the  maximum 

performance  degradation  may  not  occur  at  the  extreme  values  of  either 

T  or  f  . 

0  0 


The  parameters  to  be  specified  arc  functions  of  the  carrier 
frequency,  the  propagation  scenario,  the  link  gcorootrios,  and  the 
velocities  of  the  system  segments.  If  possible,  the  specification 
should  cover  all  of  the  possible  signal  variations  realized  by  exorcising 
the  above  factors  over  their  entire  range.  Mooting  this  specification 
with  each  link  independently  would  provide  the  utlimato  in  survivability, 
that  is,  a  system  that  can  survive  any  scenvirio  or  circumstance.  If  this 
specification  cannot  bo  met,  then  a  less  severe  specification  may  be 
possible  at  some  acceptable  loss  of  performance.  For  example,  in  nuclear 
environments,  signal  specifications  are  a  strong  function  of  the  maximum 
acceptable  outage  time.  Accepting  longer  outages  can  often  pi*ovide 
significant  specification  relief.  Additional  relief  might  be  jx)ssible 
for  systems  that  have  multiple  links  that  penetrate  the  ionosphere  at 
widely  dispersed  points.  Because  of  the  inherent  space  diversity, 
the  specification  would  not  have  to  cover  the  most  severe  case  for 
any  link,  but  some  reduced  level  of  thr'eat.  Regardless  of  the  reason, 
however,  any  reduction  in  the  specification  results  in  some  loss  of 
system  applicability. 

Methods  to  calculate  the  propagation  parameters  are  described 
in  reference  1,  excluding  absoi’ption,  which  is  covered  in  reference  2. 
Calculations  based  on  these  methods  serve  as  the  basis  for  developing 
signal  specifications.  The  first  step  in  the  specification  development 
process  is  to  choose  the  throat  environments.  For  natural  environments, 


the  threat  might  be  cither  the  equatorial  or  polar  ionosphere  during  the 
solar  maximum  depending  on  tho  location  of  the  system  links.  Nuclcoir 
threat  cnviroiunents  would  be  calculated  by  computer  codes  specifically 
designed  for  that  purpose  from  plausible  burst  yield,  altitude,  and 
location  combinations.  The  next  step  is  to  find  the  possible  link 
geometries  that  maximize  the  system  degradation  effects.  Ibis  usually 
moans  minimizing  the  angle  between  the  link  lino  of  sight  and  the  earth's 
magnetic  field  whore  the  environment  is  most  severe  and  maximizing  the 
link  path  length  through  tho  disturbed  propagation  medium.  Those  two 
criteria  can  usually  be  s.atisfiod  simultaneously  but,  if  not,  the  first 
almost  always  takes  precedence.  The  propagation  methods  in  references 
1  and  2  are  then  used  to  calculate  the  required  pro{>agation  qu;..  ‘'ties 
as  a  function  of  carrier  frequency,  area  coverage,  and  time.  The  system 
segment  velocities  are  chosen  to  provide  maximum  and  minimum  values  of 
the  signal  decorrelation  time.  Area  covoi’age  is  tyi)icany  represented 
in  contour  maps  projected  on  the  earth  of  the  proi)agation  quantities 
which  simultaneously  show  tho  area  coverage  for  multiple  values  of  those 
quantities.  In  ambient  environments,  contour  map  sets  might  be  generated 
for  different  probabilities  of  occm’rence  or  different  times.  The  nuclear 
environment  contour  sets  would  correspond  to  discrete  times  relative 
to  the  times  of  burst. 

These  contours,  whether  ambient  or  nuclear,  provide  a  data  base 
for  the  final  specification  development.  The  preferred  specification 
would  reflect  the  most  severe  conditions.  If  this  is  not  feasible,  then 
the  propagation  data  would  pemit  the  necessary  trade  offs  between  the 
threat  severity,  the  applicability  and  practicality  of  the  system. 

The  final  result  should  be  a  specification  that,  if  met,  would  provide 
adequate  system  performance  at  acceptable  cost  and  technical  risk. 


4. 


SIGNAL  SPECIFICATION  IMPLEMENTATION 


The  disturbed  signal  effects  represented  by  Equation  2  and 
the  specification  parameters  are  applied  twice  during  a  typical  system 
acquisition.  First,  the  signal  specification  is  used  during  system 
design  to  evaluate  candidate  design  solutions.  This  evaluation  is 
frequently  done  using  computer  digital  simulation  and  modeling  techniques 
which  provide  both  simulated  degraded  signals  and  detailed  descriptions 
of  the  dynamic  behavior  of  the  candidate  designs  over  all  the  specified 
signal  conditions  rapidly  and  economically. 

The  signal  specification  is  also  applied  during  testing  and 
evaluation.  At  a  minimum,  the  specification  determines  the  range  of 
effects  over  which  the  system  must  be  tested.  At  the  other  extreme, 
the  specification  may  define  not  only  the  parameter  ranges  but  also 
the  exact  form  of  the  signals.  ITiis  latter  circumstance  occurs  when 
a  system,  in  its  full  operational  configuration,  cannot  be  exercised  in 
the  ma.'iirauro  threat  environments.  For  example,  the  most  severe  natural 
scintillations  are  coincident  with  the  solar  sunspot  maximum  which 
occurs  only  once  every  eleven  years.  Systems  that  must  operate  in  nuclear 
degraded  environments  cannot  be  tested  at  all  because  of  the  1965 
Atmospheric  Te^t  "an  Treaty.  In  these  instances,  the  best  remaining 
method  is  to  degradv  the  satellite  signals  artificially  by  appropriately 
designed  link  simulators,  by  placing  a  simulator  at  one  or  more  points 
in  the  operatioru'l  system,  testing  can  be  done  limited  only  by  the 
accuracy  of  and  confidence  in  the  signal  specification. 

fclien  Uf»k  simulators  are  necessary  for  testing  and  evaluation, 
questions  arise  over  wljat  eo«istitutes  a  necessary  and  sufficient  test 
program.  For  example,  is  it  necessary  to  include  all  of  the  effects 
represented  by  Equation  2  in  every  test  simultaneously?  Is  it  necessary 
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to  test  simultunoously  every  link  in  a  system  or  even  every  link 
individually?  How  do  you  test  a  system  to  include  the  responses  of  the 
operators  as  they  react  to  degj*adation  on  system  links  from  propagation 
or  other  siraultunoous  threats?  On  a  longer  term,  is  repeated  testing 
necessary  to  insure  that  the  system  hardware,  the  operational  procedures, 
and  the  operator  res{>onscs  remain  adequate  to  maintain  the  required 
pcrfonaance?  Tlteso  and  other  similar  issues  need  consideration  for 
each  system  tested.  The  answers  to  these  questions,  broadly  known  as 
"ctrrapl iance  standards,"  should  accompany  the  signal  specification  for 
each  system  to  provide  a  completely  defined  test  and  evaluation  program. 

5 .  SUMMARY 

T!»e  preceding  sections  have  detailed  the  form  of  a  signal 

3 

structure  specification  for  application  to  satellite  (.’  systems.  The 
specification  is  intended  as  a  design  tool  as  well  a.s  a  definition  of 
the  signal  conditions  in  which  the  system  must  operate.  Adequate 
techniques  exist  to  apply  the  specification  both  in  computer  simulation 
during  design  and  by  a  link  simulator  for  test  and  evaluation. 

Appendix  A  describes  methods  to  implement  specific  realizations 
of  the  channel  impulse  response  function.  Appendix  b  provides  a  brief 
description  of  link  simulators. 
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APPENDIX  A 

THE  GENERATION  OF  DIGITAL  REALIZATIONS 
OF  THE  CHANNEL  IMPULSE  RESPONSE  FUNCTION 


The  channel  impulse  response  function  is 


h(t.T)»A(t)exp 


2Tlf^-(t) 


cr  N(t)  f 

- J  a(AO  h^U.Af) 

C  *• 


whore 


A(t)  s  absorption  and  antenna  loss 

a  carrier  frequency 

c  a  lij{ht  speed 

z(t)  a  propagation  path  length 

r  *  classical  electron  radius 
0 

N(t)  a  N, (t)*Ng(t) 

Nj It)  a  slow  trend  total  electron  content 
Ng(t)  a  yandoB  total  electron  content 


The  generation  of  digital  reaU?atior.s  of  Equation  A-l.  in 
general,  has  three  steps,  first.  A(t).  and  t(t)  are  discretited. 

Next,  samples  of  h  (t.Af)  and  Nslti  ^**0  calculated  using  Monte  Carlo 
techniques,  fijtally.  the  integral  over  Af  is  evaluated.  Most  of  the 
complexity  in  generating  these  digital  representations  lies  in  the  random 

I- 

sampling  of  h^it.Af}  and  ^^JJltj.  Most  of  this  appendix  will  coiicentrate 


on  those  sam{)ling  techniques. 


n»o  statistical  proj»erties  of  h^(t,Af)  are  defined  by 


li  (t.Af) 
s 


e 


<«*  xt 


df  h^{  f,t  )exp(  i2nAft*i2ii 


ft) 


(A-2) 


where  h^(f,T)  is  a  zero  mean  normally  distributed  random  variable  with 
its  autocovariance  defined  by 


h*(f.t) 

h.(f’,T')  .  5{f-f)MT-t*)  1'  (f.T) 

(A-Aa) 

h^if.T) 

‘'TO  »  0 

(A-Tb) 

where  <5(f)  »  Dirac  delta  function 

r-,(f.T)  «  generalized  power  spectrum 
<* 

There  are  two  distinct  cases  to  be  considered,  first,  for 
flat  fading  defined  by  f  T^l  where  f  is  the  frequency  selective  band- 
width  and  t  is  the  symbol  period,  the  generalised  power  spectrum  is 

l.Sb4T^  1(0 

I‘.(f.i)  s  p  -  -=  (A- 4) 

h*s.sr2(ff  )-r 

where  ?  is  the  signal  decorrelation  time,  because  of  the  delta 
correlation  in  delay,  h^(f.if)  is  not  a  function  of  if  and  the  problem 
is  reduced  to  generating  a  sequence  that  ts  only  a  function  of  time. 
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¥or  any  Af,  let 


h  *  Plijt.  (A-5) 

HI  H  I  *  &  I 

■>  t, 

u.  =  P  u,  J  ♦  d-p")  ^  g.  (A-6) 

where  p  a  cxp{-2.U6  At/t  ) 

'  0 

s  zevit  wean  no»TaaUy  distributed  complex  sample 
uJiTj  *  tawli  (2.l4bAt/t^) 


The  se^uence  can  be  initiated  by  choosing  an  initial  value  I'or  f;,U.,Af) 

»•  ♦ 

of  offier  unity  and  a  value  lor  U|  of  order  S*S|^  •  Appendix  b  uescribes 
a  simple  algori;S\ffl  lor  sampling  ?ero  mean  noj'mally  distributed  complex 
variables.  The  samples  Irtna  ^uation  represent  instantaneous  samples 
and  no  filtering  is  implied.  Also,  the  following  conditions  should  be 
met. 


f  /Atilb  (A*7aJ 

6 

M>iOOf  /A;  (A-Tb) 

*  6 

where  M  is  {he  miaimus  auatber  of  points  in  any  set|ue«ce  geaerated  by 
f.tjuation  A=5. 
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I 


i 


I' I 


The  secotui  case,  defined  by  i'^T^l,  is  frequency  selective  fading, 
liquation  A-2,  in  finite  difference  form,  is 


M 


h^(t.,Afj.)^AfJ_^At  r  h^(f^,tpexp(i2nfj^t.>i2t,Af^Tj)(A-8) 


where 


Af  »  l/(2MAt) 

*'  ■  'i-'i-i 

f  =  n.A  f ' 
n 

V  iAt,  iil<M 
|Afj.l  <  1/(2At) 


From  Equations  A- 5a  and  A-Sb, 
be  written  as 


the  statistics  of 


?  fi*At/2 


b*(f  .f  i)h  (f  .'f  . )  s  3  3. 

s  n  j  s'  a*  i  lan  *nAi 


^  t  .At/'2 


/ 


f„*AfV2 

df 


Mf.f) 


tA-ga) 
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h. 


't 

t' 


i 


.< 


whero,  for  f  T<1, 
0 


f  =  t-^dK-J)’* 

Cj  =  delay  {wroweter  (^0.2S) 

Iquation  A- 10  and  related  quantities  are  evaluated  in  Appendix  C. 
tlsnerating  the  h  (t.,Af.)  consists  of  evaluating  the  h.(f’  .^•)  using 
liquations  A-9,  A-IO,  and  Appendix  B  and  then  evaluating  Equation  A-8. 
For  adequate  statistical  sampling  and  numerical  resolution,  it  is 
necessary  that 


T  /AtS'lO 

0 

(A^Ua) 

/At 

h 

(A-Uh) 

At<T/S 

(A. tic) 

AlAf)<i/aax(l.l/f  .lOT) 

0 

(A-ild) 

h  f[<0.75 

G 

(A. tie) 

-0.2S<2«f  t<S.4S 

tA-tlf) 

where  ^ .  Condition  A-llb  is  a  bate  miniwum.  I.arger  M  or 

multiple  sequences  are  advisable,  liquations  A-lle  and  A-llf  prescribe 
the  ranges  of  t  and  f  in  liquation  A- 10  tltat  include  9,^  percent  :>f  the 
total  energy.  The  coefficients,  ‘'of  necessary  outside 

those  ranges,  liquation  A-llg  and  the  T  dct>endcnt  tern  in  liquation  A-ild 
reflect  an  estimate  of  the  required  delay  range  to  handle  dlsiK?rsion. 

TTie  integral  in  liquation  A-9  can  he  evaluated  to  one  percent  with 
trapezoidal  integration  with  frequency  and  delay  increments  less  than  or 
equal  to  0.03/t^  and  U.05/f  ,  respectively. 

The  representation  used  for  the  random  portion  of  the  totai 
electron  content  depends  on  the  application.  The  first  use  of  in 

liquation  A- 1  is  in  the  carrier  phase  where  it  can  degrade  phase  or 
frequency  acquisition  and  tracking,  for  pluse  effects 


M* 

)=Af**  I  g(f ’)expCi2nf’t .  )  (A*12) 

ns-M*  ’  "  ♦ 

where  Af”  =  i/(2M'At‘) 

f  *  =  ni  f  “ 

ft 

t .  s  iA{ *. I i|«M- 

g(fj|)  is  a  sero  -sean  floi^^^ily  distributed  variabie  whose  resainiag 
statistics  are  defiaed  by 


f>At'’72 


r*i» / 
')  •'f*.Af*72 


where  g(f*  )  »  ) 

-n  **  '  n 

“  "Wr'/ 

fg  »  Ra/leigh  tVe«|ueney 

{aiUal  estimatesi  fov  At*  aji^  M*  are 

At*<l/(4ty 

M*>$U* 

for  ti'e^ueney  traekiag.  Ec|uatieo  A-U  eaaaat  he  used,  tas^ead.  let 


■ 

(A-I5a) 

-  0 

(A-lSb) 

(A-Ub) 


(A  =  ||) 


The  second  use  of  N(t)  in  Equation  A-1  is  in  the  dispersive 
phase.  Nj^(t)  can  usually  be  omitted  because  the  dispersive  effects  are 
primarily  determined  by  the  magnitude  of  N(t)  and  not  the  time  dependence. 

The  last  use  of  N^(f)  determines  the  group  delay  and  group  delay 
rate  which  can  degrade  time  synchronization  if  sufficiently  large. 
Equations  A- 15  and  A- 16  should  be  used  with  fj^  replaced  by  f.^.. 


i 

I 

I 

] 

a 

I 

I 


fT=f„/  '' 


(A- 17)  I 


If  simultaneous  handling  of  the  total  electron  content  effects  is 
necessary,  then  Equations  A-15  and  A-16  ai’e  necessary.  Tlie  substitution 
of  f.p  should  not  be  used.  Unfortunately,  the  number  of  time  increments 
necessary  to  span  the  time  range  can  be  very  large.  Exercising  the  range 
of  time  may  not  be  necessary  in  all  cases.  Some  experimentation  is  often 
indicated.  Another  problem  with  Equation  A-15  is  that  the  sura  cannot  be 
evaluated  with  fast  fourier  transforms.  Finally,  multiple  sequences  of 
Nj^(t^)  should  be  used  to  test  the  adequacy  of  the  sampling  for  either 
Equation  A- 12  or  A-15. 


The  entire  channel  impulse  response  function  is 


2:Tf  2(t  )  cr  N(t  )  I  ^ 

h(t.  ,T  .)=A(t.)exp  i - - — i--i - j -  A(Af)  I  h^(t..Af|^) 


crN(t.)AfJ  .  /  2(tjJ  cr^N(t,) 

exp  -i - -  i^irAfjT.-  ^ 


(A- 18) 


m 


where 


L=l/(2ATA(Af)) 


cr  N(t.) 

0  1 


<LAr 


In  principle,  one  only  needs  to  insert  the  proper  quantities 
discretized  to  the  smallest  increment  in  time,  delay,  or  Af^,  and  suia 
over  At'j..  In  practice,  however,  the  wide  range  of  the  delay  and  time 
necessary  to  adequately  sample  the  various  functions  make  this  impractical. 

It  is  usually  necessary  to  handle  some  of  the  terns  indejiendently , 
particularly  the  total  electron  content  and  path  length  determined  effects. 

The  channel  impulse  function  is  thereby  greatly  simplified. 

h(t.,tj)=A(tph^(t^,t.)  IA.19) 

where  for  f  Tsl 

M 

£  h  (f  ,r  )exp(i2flf  t.) 

or  for  f«T5“l 

(5 

h^(ti.'f,).h^(t..Af).  T.=() 
sO  . 

A  fortran  program  for  generating  h  (t-.T-)  can  be  found  in  Appendix  E. 

Equation  A-19  is  implementable  in  a  link  simulator  as  described  in  Appendix  0. 

An  important  consideration  for  implementing  Equation  A-18  or  A-19 
is  the  dynamic  range  requii*ed  to  resolve  the  JiapUtude  fluctuations,  if 
present.  The  amplitude  of  the  channel  impulse  response  function  at  each 
discrete  delay  is  a  Rayleigh  distributed  variable,  let  the  amplitude  of 
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be  R,j.  The  distribution  function  for  R..  is 


2R.  .  /  7  oX 

P{K.j)=-^  oxp^-R.j/oJj 


(A-20) 


where 


'■/"I* 


t.Mt/2 


liquation  A-20  is  easily  invertable  to  form  probability  statements  about 


U. ^ .  thus, 


I -exp 


■(fl 


» probability  that 


s  probability  that 


We  require  that  the  amplitude  not  exceed  the  dynamic  range  99.8  percent 
of  the  time,  the  upper  and  lower  bounds  on  the  amplitude  are, 
respectively 


®  j^-ln(0.00l)j  5  2. 65 

^  '  *  5,16; 


16x10  ■ 
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where  the  protnibil  ity  that  the  amplitiKle  exceeds  the  dynamic  range,  0.002, 
is  equally  split  t)ctween  the  limits.  In  decibels 


30dB<20 


Thus.  .^8.4  dli  of  dynamic  range  is  required.  If  this  is  not  feasible,  then 
the  loss  of  accuracy  should  be  at  the  largest  amplitudes,  l*)»is  is 
because  system  performance  is  much  moi*e  sensitive  to  amplitude  fades  than 
to  amplitude  enliancements . 

The  resolution  in  the  amplitude  and  phase  is  a  functioJi  of  the 
resolution  of  the  real  and  imaginary  cotqionents  of  b(tpt^).  The  {Joorest 
resolution  occurs  in  deep  fades  at  the  smallest  amplitude.  The  maximum 
increment  is  defined  by  requiring  that  successive  values  of  either 
comi)onent  of  htt..tj}  not  vary  by  more  than  ten  percent.  This  gives 
adequate  accuracy  on  the  amplitude  and  resolves  the  phase  to  within  six 
degrees.  If  the  component  increment  is  fixed,  then  the  smallest  value  must 
be  used.  The  fixed  increment  for  either  component  can  be  expressed  as 
one  tenth  of  Variable  increments  are  more  desirable 

because  they  permit  more  efficient  storage  of  cataponent  data.  They  also 
pei-mit  the  use  of  logarithmic  digital  to  analog  converters  in  link 
simulators  as  described  in  Appendix  0. 

The  total  electron  content  effects  neglected  in  liquation  A- 19 
can  be  bounded  to  determine  necessary  but  not  sufficient  conditions  for 
successful  sysftr-ffl  operation.  The  total  electron  content  and  its 
derivatives  are  gaussian  distributed  random  variables  with  the  following 
means  and  standard  deviations. 
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d”N(t) 

dt® 


s 


(A-21a) 


0 


Nm  ^ 


(A-2n>) 


These  statistical  moments,  when  converted  to  moments  of  phase,  doppler, 
dopplev  rate,  jerk,  etc.,  allow  perfomance  dejjradutlon  estimates.  These 
statistics  also  determine  the  amount  of  time  that  any  of  the  derivatives 
exceed  some  value.  As  such,  they  provide  the  means  for  specifying  various 
system  parameters  consistent  with  some  perfoiwnce  requirement. 


I'll©  total  electrcHj  content  drives  the  dispersive  distortion  of 
the  signal  waveform.  As  a  measure  of  this  distortion,  we  use  the  expec¬ 
tation  value  of  the  absolute  value  of  the  delay  coordinate. 


«l3 

f 

oa 

w 

cr  NitlAf" 

J  dr  !  r  [ 

dCafjexp 

- -i2«AfT 

r 

c  J 

CA.22a) 


B 


cr  Ni(t) 

6 


“ITT 

c 


(A.2|b) 


Since  the  statistical  moments  of  the  measure  of  distortion  are  proportional 
to  the  m<»ents  of  the  total  electron  content,  performance  criteria  can 
be  developed  as  a  function  of  modulation  and  other  system  details,  the 
criteria  should  guarantee  tt»at  the  ratio  of  i?!  to  T  he  small. 


The  finite  difference  equations  have  been  written  in  the  lowest 
order.  This  was  done  for  simplicity  and  to  be  consistent  with  fast 
fourier  transform  algorithm:;.  Those  algorithms  can  evaluate  liquations 
A-8,  A-12,  A-18,  and  A-19  much  more  efficiently  than  conventional  methods. 
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APPENDIX  B 

A  RANDOM  NUMBER  GENERATOR  FOR 
NORMALLY  DISTRIBUTED  COMPLEX  NUMBERS 


Let  g  be  a  nomally  distributed  cofli{>lox  number  defined  by 


— 2 

8*8  »  (U-Ia) 

gg  =  U  (B-lb) 

8  *0  (B-lc) 

i’(lgi)  •»  4^*exp(-!g|  “/a”)  (B-ld) 

«• 

0 


where  PClgl)  is  the  pi‘obabtHty  distribution  for  the  amplitude  of  g.  Let 
be  a  unlfomly  distributed  random  variable  on  the  iatet^val.  t><R^<l. 

A  sample  of  g  is  ealeulated  by 


g  ’ 


a|^'ln(Rj)j  ^eJtp(i2»R^^j) 


(B-’) 


APPENDIX  C 

EVALUATION  OF  THE  GENERALIZED,  DELAY. 
AND  FREQUENCY  POWER  SPECTRUHS  FOR  f  T<I 


The  gonoralized  power  spectrum  is 


where  f 

L‘j  *  delay  paraaeter  (»o.JS) 

?  frenuency  selective  bandwidth 
fg  '  signal  decarrelat ian  tlae 

A  nt  ta  Kquatiaa  C-l.  gaad  ta  ane  percent,  is 


3J 


c.(z)  w  ^  0.II21 


2  2  2 


fJ(2)  »  1 .8l3cxp^-0.(,752.U.7292'^-0. 


1092  ^0.031 


-  O^.jTj,  O.OOlA 

\z\\  ‘  2^  ^“7* 


tUi*  U©iay  power  spewtrusi  is  Jofinoil 


Ut 

at* 


f-^Cr)  5  ni’ 


« 


t  ?  {i*0.4?&i7f 


r^(t)  a  Tjf’t'xpi 


(2nf'T) 
2 


2C 


{■•(t)  .  2S() 


(C-Sa) 


r,(T)  *>  -2pf’t^|^2-{’(t)tjx|>{-r^)^  ,  2^0 

t’(t)  a^(0.7478SS6t-O.U9Sii79S)t  ♦  0.5480242^ 

Kquatlon  C-5  is  accurate  to  about  five  figures.  The  frequency  {Hwer 
spectrum  is  defined  as 


(C-5b) 


Tjlf)  « Jdtr^lf.t) 
r^(f)  *  p'f^exp|^'(»t^fj'‘ j 


(C-b) 


(C.7) 


The  normal  isatiea  integrals  for  the  abov^e  spectruias  are 


u.  a. 


H 

A 

/ 


df  -  I 


(C*8) 


it)  »  I 


iC.9) 


dfrjf)  5  l 


fC-iO) 


ss 


Frequently,  integrals  over  the  various  power  spectrums  are 
necessary  as,  for  example,  in  L-quation  A-d.  Trapezoidal  integration  with 
the  following  integration  increments  is  accurate  to  approximately  one 
{lercent . 

dfj  i  0.03/t^^  (C-lla) 

ilTj  i  0.03/f^  (C-llh) 


APPENDIX  0 

PROPAGATION  EFFECTS  LINK  SIMULATORS 


P i ^ure  n- 1  shows  a  gottorat  iin{>teoicntatian  oF  a  link  simulator 
which  is  the  hardware  aualoii  of  the  chamtel  impulse  res{M>nse  fuiKrtiun 
defines!  by 


«(t)  »|  dth{t,T)S(C-T) 


{P-D 


^1  ! 


wRere  i»lE)  »  transmitted  linput)  signal 

R(t)  55  received  (output)  signal 
htt.T)  »  chattel,  impulse  response  function 


the  link  simulator  represents  a  discretised  version  of  P-l  where 


yt)  . 


CB.:a) 


q^Ct)  s  P/.A|laag|^htt. j| 


P/a{  I  3  digital  to  aaalof  couvifrstoa 

Xlgorithsi  to  calculate  rvalues  of  tne  chanfwjl  liopulse  response  fuaction 
were  detailed  is  Appeedi*  A.  tfee  sailing  t‘e^i?eeeais  iu  Appeedis  A 


F-.gure  D-1 .  Propagation  Effects  Simulator 


are  sufficient  to  support  the  digital  to  analog  operation  without  any 
additional  filtering  to  suppress  aliasing. 

Figure D-2  shows  a  particularly  simple  link  simulator  for  flat 
fading  and  when  all  other  effects  are  ignored  or  simulated  separately. 

The  delay  is  used  to  get  an  independent  noise  source  for  the  1 
channel  without  an  extra  noise  generator.  The  filters  are  single  pole 
providing  in  cascade  the  power  spectrum  represented  by  Equation  6.  For 
a  RC  filter,  the  time  constant  is 


RC  =  T  /2.146  (D-3) 

0 


The  input/output  can  be  at  any  frequency  level. 


APPENDIX  E 

CIRF,  A  FORTRAN  PROGRAM  FOR  GENERATING  SAMPLE  DIGITAL 
REPRESENTATIONS  OF  THE  CHANNEL  IMPULSE  RESPONSE  FUNCTION 


This  appendix  contains  a  fortran  program  to  sample  digital 
representations  of  the  channel  impulse  response  function  in  Equation  A-19 
less  the  absorption  and  antenna  loss  terms.  The  coding  is  consistent  with 
Fortran  IV  with  the  exception  of  the  dimension  declaration  in  the  fast 
fourier  transform  subroutines  (see  comment  in  SUBROUTINE  SETD  list). 

The  methods  used  to  generate  the  frequency  selective  sequences 
are  also  used  for  the  flat  fade  sequences  rather  than  Equations  A-5  and 
A-6.  The  flat  fade  sequence  generation  is  thus  slower,  but  the  sequences 
generated  are  continuous  from  the  last  sequence  term  back  around  to  the 
first  term,  a  very  useful  property  for  hardware  testing  applications. 

The  flat  fade  sequences  are  nevertheless  consistent  with  Equation  A-4. 

The  program  tests  the  input  for  compliance  with  Equations  A-Ua, 
A-llb,  A-lle,  and  A-llf.  Tlie  recommended  integration  increments  are  also 
enforced.  The  use  of  floating  point  variables  for  the  channel  impulse 
response  function  output  guarantees  compliance  with  the  stated  phase  and 
amplitude  resolution  requirements  on  all  known  computer  systems  that 
support  Fortran  IV.  The  user  is  responsible  for  satisfying  Equation  A-Uc. 

The  code  provides  time  sequences  of  the  impulse  response  function. 
The  first  sequence  centers  on  a  delay  of  .0.04/l'0*PD/2  where  FO  is  the 
frequency  selective  bandwidth  and  00  is  the  delay  increment.  The  main 
program  loop  over  delay  ends  at  statement  label  70.  Just  prior  to  that 
statement,  the  time  sequence  for  the  current  delay  is  in  the  H  array. 
H(I>I-1)  and  represent  the  in-phase  (real)  and  quadrature  (Imaginary) 

impulse  response  function  co^>onents,  respectively.  The  user  is 


responsible  for  providing  the  extra  coding  for  disposing  of  the  sequence 
data.  The  data  must  be  moved  out  of  the  array  before  the  loop  returns  to 
calculate  the  next  sequence. 

Any  remaining  questions  should  be  referred  to 
Director 

Defense  Nuclear  Agency 
RAAE/Major  Leon  A.  Wittwcr 
Washington,  DC  20305 

Autovon  221-7028 
Commercial  202/325-7028 

The  program  was  programmed  and  debugged  on  a  Hoathkit  H89 
Computer  system. 


PROORAM  CtRF 
C 

C  THIS  PROGRAM  CALCULATES  RANDOM  SAMPLES  OF  THE  SCINTILLATION  CHANNEL 
C  IMPULSE  RESPONSE  FUNCTION  AS  DEFINED  IN  "A  TRANS-IONOSPHERIC  SIGNAL 
C  SPECIFICATION  FOR  SATELLITE  C**3  APPLICATIONS."  (ENCLOSING  DOCUMENT) 

C 

C  THE  PROGRAM  INPUTS  ARE 
C  TAUO.  SIGNAL  DECORRELATION  TIME (SEC) 

C  FO,  FREQUENCY  SELECTIVE  BANDWIDTH(H2> 

C  OT,  TIME  INCREMENT (SEC) 

C  (THE  TIME  INCREMENT  MUST  BE  LESS  THAN  OR  EQUAL  TO  TAUO  DIVIDED 

C  BY  TEN) 

C  DD.  DELAY  INCREMENT (SEC) 

C  (THE  DELAY  INCREMENT  SHOULD  BE  LESS  THAN  OF  EQUAL  TO  THE 

C  RELEVANT  SYSTEM  SYMBOL  DIVIDED  BY  FIVE.  IF  DD  IS  GREATER  THAN 

C  OR  EQUAL  TO  SIX  TENTHS  DIVIDED  BY  FO.  THEN  ND(SEE  BELOW)  IS 

C  SET  TO  ONE  AND  THE  FLAT  FADE  SPECTRUM  IS  USED) 

C  NT.  NUMBER  OF  DISCRETE  TIMES 

C  (NT  MUST  BE  A  POWER  OF  TWO,  NT  MUST  BE  SUFFICIENTLY  LARGE  SUCH 

C  THAT  NT  TIMES  DT  IS  GREATER  THAN  OR  EQUAL  TO  ONE  HUNDRED  TIMES 

C  TAUO) 

C  NO.  NUMBER  OF  DISCRETE  DELAYS 

C  (NO  MUST  BE  sufficiently  LARGE  SUCH  THAT  ND  TIMES  DO  IS  GREATER 

C  THAN  OR  EQUAL  TO  SIX  TENTHS  DIVIDED  BY  FO) 

C  SEED.  SEED  FOR  RANDOM  NUMBER  GENERATOR 
C  (SEED  MUST  BE  AN  ODD  REAL  NUMBER) 

C  REQUIRED  FUNCTIONS/SUBROUTINES 
C  FLOAT  (FIXED  TO  FLOATING  POINT  CONVERSION) 

C  SETD  (INITIALIZES  FAST  FOURIER  TRANSFORM  TABLE) 

C  IFIX  (FLOATING  TO  FIXED  PCUNT  CONVERSION) 

C  DINT  (DELAY  INTEGRATOR) 

C  DFAM  (DF*(STANDARD  ftiVIATION  OF  FOLIRIER  COEFFICIENT)) 

C  RANCGF  (RANDOM  DFiKFOURlER  CCEFFIC lENT )  GENERATOR) 

C  PFT  AND  PDO  (FAST  FOURIER  TRANSFCiRM) 

C  PO(«tER  (SIGNAL  FLOWER  DENSITY  AS  A  FUNCTION  OF  DELAY) 

C  TCOR  (COMPLEX  SIGNAL  TIME  OECORRELATION  FUNCTICW) 

C 

DOUBLE  PRECISION  SEED 
C 

C  THE  rOLLOUING  DIMENSION  LIMITS  AR£  MINIMUM  VALUES 
C  DIMENSION  H(2«NT>.D(NT> 

c 

EUMENSieiN  HI  40'5'R> .  0120*9) 

COMMON  /DAT/NFP.OFP.NPP.DOP.ND.M* 

COMwriN  /rsD/  e  I .  T auo  .  po 
eOHM>?N  /RAM/SEED 
COMMON  /tiFM/DFMS 

c 

C  SET  delay  parameter 

c 

DATA  ei/.2<4/ 

C 


C  READ  AND  TEST  DATA 
C 

C  INPUT  IS  ON  LOGICAL  UNIT  6 
C  OUTPUT  IS  ON  LOGICAL  UNIT  2 
C  OUTPUTICRTJ  IS  ON  LOGICAL  UNIT  1 
C 

CALL  0PEN<<..'SY1J  INPUT.DAT  ') 

CALL  OPEN <2. 'SYltOUTPUT.DAT  ') 

READ! 6 . 1 000  >  TAUO . PO . DT , DD. NT . NO. SEED 

1000  F0RMAT(4£10.2.2I5.D10.2) 

WR I TE < 2 . I 00 1 ) T AUO . FO . DT . DD . NT . ND . SEED 

1001  F0RMAT<IX.1P.27HSI0NAL  DECORRELATION  TIME-  .ElO.2/  IX.SIMFREOUENCY 
1  SELECTIVE  BANDMIDTH-  ,E10.2/1X. 16HTIME  INCREMENT-  .E10.2/1X. 17MDE 
ILAY  INCREMENT-  .E10.2/1X. 17HNUMBER  OF  TIMES-  . IS/IX. leMNUMBER  OP  0 
lELAYS-  .I5/1X.6HSEED-  .020.12) 

URlTE(l.tOlO) 

1010  FORMATItX.SHRUNNINO  > 

IF(DT.LE.TAU0/10. )00  TO  10 
WRITEiS* 1002) 

1002  FORMAT! IX. 33HERR0R-T1ME  INCREMENT  IS  TOO  LAROE  > 

STOP 

10  AMP-FLOAT < NT) 

NP-0 

20  AMP-AMP/2. 

I 

IF(AMP.0T.1.)00  TO  20 
IFtAMP.EO. I. ICO  10  30 
WRlTt:<2.1003> 

1003  FORMAT!  IX.  3>>HERR0R-NT  IS  NOT  A  POUER  OP  TUO  ) 

STOP 

30  IF!FLOAT!NT>*OT.OE.100.*TAUO)00  TO  «0 
URITt!2.1i>04) 

1004  FORMAT! IX. 3SHERR0R'-TIME  RANGE  LESS  THAN  100  TAUO  ) 

STOP 

40  IF!DD*FLOAT<ND).C«E..«./FO)00  TO  50 

URITS’!?.  IO<S^) 

1005  FORMAT! tX.34HERR0R-DELAV  RANGE  LESS  THAN  O.A/FO  ) 

STOP 

50  IF!PP.LT..4/FO)00  to  40 

ND«| 

URITE!2.  l€iOA) 

1004  FC«MATt)X.39HFLAT  FADE  LIHIT-ONLV  ONE  DELAY  K£OUIK£D  > 

40  CONTINUE 
C 

e  initialiie  fast  fourier  transform  tadle 

c 

CALL  £ETD!O.NP) 

L-NT*NT-t 

e 

t  calculate  FRCOUENCV  increment  and  INVCORATION  VAR|ADl£S 

c 

e  pPeFREOUEKCV  increment 
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Of»l./(FLOAT«NT)*OT) 

C  NFP«NLiri6EFl  OF  FREOUFNCV  SUBINCREHEMTS 
NFP«lFI*(33.«TAl.iO»0F»^» 

C  tlFP-FREOHENCY  SUP  INCREMENT 
DFP-OF /FLOAT «NFP> 

C  NOP-NUMPER  OF  OELAV  SUP INCREMENTS 
NDP-IF I X ( 33. *00»FO>*t 
C  OOP-DELAV  SUB INCREMENT 
00P-DD/FL0AT(NDP) 

C 

C  CALCULATE  TIME  SEQUENCES  FOR  EACH  DELAY 

C 

C 

C  DSTART  IS  THE  BEOINNINO  OF  THE  t€LAY  WINDOW.  TWO  PER  CENT  OF  THE 
C  signal  ENERY  arrives  before  DSTART  AND  THUS  IS  NEGLECTED  WHEN  DSTART 
C  IS  INITIALUEO  TO  -.04/F0. 

C 

DSTART*-. 04/F0 
C 

C  INITIALIIE  VfARIABLES  FOR  STATISTICAL  TESTING.  THE  AVERAGE  SIC-NAL 
C  POWER  AT  EACH  DELAY (THE  AVERAGE  VARIANCE  OF  THE  FOURIER  COEFFICIENTS 
C  times  00*«2)  and  the  TIME  DECORRELATION  PROPERTIES  OVER  ALL  SEQUENCES 
C  NEAR  TALtO  WILL  BE  TESTED. 

C 

ICOB»IFIX(TAUO/DT) 

icoR»ieoR*icoft*i 

CORRsO. 

CORl»0. 

DO  70  1*1. ND 
T»OSTART*.S«DD 
WRITFJ  J.tOUIT 

1012  FORMAT <  I IP.  JIMGENERATINO  SEQUENCE  FOR  DELAY*  .CI0.2> 

WH. 

F?TART».S*eP 

AMP»l. 

C 

C  ZERO  DC  FREQUENCV  COMPONENT  TO  MAINTAIN  ZERO  MEAN  SEQUENCE 

e 

Ht|t«^. 

N(2I*0. 

c 

c  DFHSBlNlTtAL  (‘ELAY  INTEGRATION  POINT  FQR  DPAH 
C 

DFHS^.f  ♦DINT  ( EiSTART  .FSTAST  » 

DO  60  ytejI.NT.S 
IFtAMP.LE.O.  fO0  TO  6?. 

AMPsDFAMtDSfART.FSTAftT » 

CALL  RANCOP(H(Jt.M(J*|l.AMP» 
call  PAKCertHCT  ).Ht>.*tt.AMPl 
FSTAflT?^':TAftT«eP 
oa  TO  eo 
isS  N(UI«0. 
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H<K)«0. 

H(K*l )»0. 

SO  K-X-2 

AMP-PFflM  <  DiTART , FSTART ) 

CALL  RANCOF(H<K),H<K*l).AMP) 

C 

C  use  FAST  FOURIER  TRANSFORM  TO  GENERATE  FINAL  SEQUENCE 
C 

CALL  PFT<H.D> 

CALL  FDD<H) 

C 

C  TEST  SIGNAL  POWER  IN  SEQUENCE  TERMS  AND  ACCUMULATE  DATA  FOR  TIME 
C  CORRELATION  TEST.  USE  TRAPEZOIDAL  iNTEORATlON. 

C 

TVAR«0. 

DO  90  U-IC0R.L.2 
K»U-ICOR^l 

CORR=CORR»H<K>*H(U)*H(K>I>*H(U*1) 

CORI»CORI'»H(K*n*H«J)-M(K)*H<U*t> 

00  TVARaTVAR*H< J>*H(U>>H(U*I)*H(J*l) 

AMP>«2 . 0*  T  VAR*DC**  •  2/  FLOAT  ( L-  IC0R*2 ) 

PSUM« . 5*P0WER  <  DSTART ) 

T-OSTART 
DO  too  Uat>NOP 
T»T*DDP 
SM^POUERIT) 

100  PS.UM«F";.UN*$M 

PSUM» ( PSUM* . ?*SM ) •DDP 

T»DSTART*.?*PO 

UR  I T  E  <  2 .  I  OOr? )  T .  AMP .  PSUM 

lOOC  FORMAT* / 1  X.  tP.7MPELAV=  .£10.2/1  X . 24HAVERAC.E  SEOtjCNCe  PQUERs  . 

lEt0.2/U.  i^-sHTHEORETICAL  SE'SL^NCE  POUER*  .£10.2) 

C 

C  THE  H  ARRAY  NOW  cWTAlNS  A  SAMPLE  DIGITAL  REPRESENT  AT  ION  OF  THE 
C  CMANfJFL  IMPULSE  Rf:i-PON-;E  FUNCTION  FOR  A  DELAY  OF  £i'START*.S*C't'.  THE 
e  IN-PtiASE  AND  OUApRATURg  W^iLTlPLlERSCTHg  REAL  ANp  IMAGINARY  PARTS  OF 
C  THE  CRANNgL  IMP'JLSE  RESPONSE  FUNCTION.  REiPECTlUELV)  ARE  DEFINED  DV 

e 

C  IN-PHASC  eOMP*5NENT  MyLTIPLlFR  s  Htl*I^t) 

e  QUAt>f<ATgRE  COMPANY  MULTIPLIER  *  M«I*I> 

e 

e  WHERE  I  INi>E*F;S  SCtCCESSIVE  VALUES  IN  TIME.  THE  U'JER  MUST  rCCUC  ON 
C  TF*E  plSPOSlTtON  OF  TmE  L<AIA  IN  TMg  H  AF:RAV  egFORg  PROCEEDthO  ON  TO 
e  THE  HE»T  tCLAY.  USUALLY.  THE  DATA  IS  STOREf)  AND  LATER  USED  IN  AN 
e  IHTEORAL  O'/fR  DELAY.  IN  These  iNTEORATJijN'?.  DD  MUST  £*£  USED 
C  eiPLieiTLV  AS  THi  INTEORATtVN  INCREMENT. 

c 

VO  DSTARTsDSTART^tiD 

e 

e  TEST  TIME  eORRLLATtUN  PROPERTIES 

e 
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o  <ti  rt  .n»  fl 


COfiR-2.  rt*CORR*tiD*  •2/FUOAT  <  1.  - !  C0R»2  > 
CORl«2.0*CCiRt*DD**2/FL0AT(L-:C0R*2> 

T». 5«FL0AT  t  tCuR- 1 )*DT 
AnP«TCCIR(T) 

WRlTe«2. tOO«>T,CORR.CORI.AMP 

J009  F0RMAT</1X.  1P.24.HTIME  D!?.FuACE«ENT  TESTED-  .EtO.2/  lX.26MCAl,ait.ATE 
10  DECORRELATION-  .2E10.2/lX.25HTH£0ft£T:CAC  CALCULATION-  .E10.2> 
WRITE(2. 1007) 

1007  FORMAT! U.IOHENO  OF  RUN  » 

END 

FUNCTION  OFAM<DSTART.FSTART) 

C 

C  THIS  FUNCTION  CALCULATES  DF  TIMES  THE  STANDARD  DEVIATION  OF  THE  RANOC'M 
C  FOURIER  CORFFICIENT  CORREiPONDINO  TO  A  OELAV  OF  DSTART*. S-DD  AND  A 
C  FREQUENCY  OF  FSTART*.5*0v“.  THE  VARIANCE  OF  THIS  FOURIER  COEFFICIENT 
C  IS  THE  GENERAL  I  ZED  POWER  SPECTRAL  DENSITY  InTEORATED  OVER  DELAY  FR.jM 
C  OSTART  TO  DSTART+DD  AND  OVER  FREOUENCY  FROM  PSTART  TO  FSTART+DF, 

C  INTEGRATION  IS  BY  THE  TRAPEZOIDAL  RULE. 

C 

C  ARGUMENTS 

C  DSTART- INITIAL  DELAY! SEC ) 

C  FSTART- INITIAL  FREOUENCVIHZ ) 

C  REQUIRED  COMMON  VARIA&LES 
C  NFP-NUMPER  OF  FREOUENCY  SUBINCRE1«NTS 
C  DFP-FREOUENCV  SUDlNCREMENTiMZ) 

C  ODP-DELAY  SUB INCREMENT (SEC) 

C  DD-DELAV  INCREMENT 
C  DFMS-0. 5-DlNT  <  DSTART . FSTART  > 

c  Functions  required 

C  DINT  !DELAY  INTEGRATOR) 

C  SORT  !  float  INC*  POINT  SCsUARE  RQiOt) 

c 

COMMON  /DAT/NFF‘.DFP.NEtP.OCP.ND.DO 

COMMON  /t'PM/DFMS 

F-FSTART 

Df  AM=DFM-3 

DO  10  Isl.NfP 

F*F*DFP 

SM^DINTCDSTART.P) 

10  pFAM=fiFaM*S« 

DFMt=.S*SM 

fiF  AM=iOftT !  DFP-DDP*!  DFAM-DFMS  n /DD 

return 

END 

FUNCTION  DINT ( DSTART. R) 


THIS  FLiNeTtON  INTE0P6TFS  TmE  OfNER&LlJgD  EfeoH  fe  DELAY 

OF  DSTASf  to  DSTART.DD  at  A  FR£'X€NCV  OF  F  AN©  DIVIDES  fwE  RESULT  &V 

eop. 


AROtiMfNTS 

DifART^INITIAL 


uuu 


c  r>rREauENCV(HZ) 

C  RCOtURED  COMMON  VARJAfiCES 
C  NO-FLAT  FADE  FLAO 
C  NO-I.FLAT  FAOtNO 
C  NO.OT. l.FREOLCNCY  SCLECTtVE  FADING 
C  NOP-NLIMBeR  OF  DELAY  ,?U8 INCREMENTS 
C  DDP-OELAY  SUBINCREMENT 
C  FLWCTIONS  REQUIRED 
C  PS03  {GENERAL  I  ZED  POUER  SPCCTRUH) 

c  RSOFF  {Flat  fade  spectrum) 
c 

COMMON  /oat/nfp.off»nop.oor.no.oo 

TEST  FOR  flat  FADING 

IFtND.EO.nOO  TO  tOO 
O-OSTART 

DINT-.5*PS03<F.D) 

DO  to  I«t.NOP 
0-0*DDP 
SM«PSP3{F.0I 
10  0INT-01NT*SM 

OINT»OINT-.S*S« 

RETURN 

too  OINT-PSOFF{F)/OOP 
RETURN 
END 

function  PSDTIF.D) 

c 

C  TM!S  FUNCTION  1$  THE  GENCRALIZEO  POUER  SPECTRUM  FOR  FR£)^JCNCV 
C  SELECTIVE  CASES. 

C 

C  AROUWEN'S 
C  F*FREO»^NCY«HZ) 

€  0=I>C;LAV<«ge) 
e  REOfJIRE^  C0MM»3N  VARIAfiLES 
C  et-DELAV  PASAMETER 
e  TA»J«>=SlONftL  DECOSRELATION  TIME 
c  FQ^FREOUENCV  1  elective  |iAn£<Ui{OTH 

e  fvni;tion‘5  rci^^jiked 

C  S((ST  < float iN'i  point  SClMAftE  ROOT) 

C  ABS  trLOafSfiO  POINT  AS'-OLOtS  VALL^E » 

C  EDP  {FLOATING  FOINT  NATURAL  E«PONENTtAL) 

c 

COMMON  /PSO/Cl.TAUO.FO 
FPRtME^O*'-.i?RT( 

TPFPO^fc.  ?:S5*FP{ifME*0 

{Tt.  tAS^TALfi^Ft-^S-TPEPOiA:!) 
lF{A|tS{T).LT.t.»«>G  TO  Id 

IFfZ.LT.O. >60  TO  20 


o  o  o  o  o  o 


00  TO  50 

20  Oi«l.772*<l.*«(.OOt4«X-.l78)***.29»*X)/SOI»T(-Z) 

GO  TO  SO 

10  02»<  <  < .031*2“. 109)*Z-.729)*2-.67S)*2 

0Z-I.3I3*EXP(G2) 

50  PS03*2.981*FPftlMe*TAu0*£XP<.5*Cl**2-TPFP0>*02/S0ftT«.l » 

RETURN 

END 

EUNCTtON  PSOEFiF) 

C 

C  this  function  is  the  flat  fade  POWER  SPECTRUM.  THIS  FUNCTION  IS  NOT 
C  THE  OENERAUtZEO  P*3WER  SP£CTRUMtFw:TtCiM  PS03>  INTE0RA7E0  OVER  ALL 
C  OELAV.  A  CLOSE  APPROXIMATION  HAS  SEEN  CHOSEN  INSTEAD  BECAUSE  OF  ITS 
e  utility  in  SOFTWARE  AND  mARDNARE  APPLICATIONS.  WHITE  NOISE  FILTERED 
C  BY  TWO  SINGLE  POLE  FILTERS  REPRODUCES  THIS  SPECTRUM. 

C 

C  ARO».'MENTS 
C  F-FRE'XCNCYiMZJ 
C  REOUIRED  COMMON  VARIABLES 
C  TAOO*SIONAL  decorrelation  TIME 

c 

CCtMMON  /PSO/Ct  .TAijO.FO 

PSOFF* 1 . 0A4*TAU0/ ( t . *0. 572* ( F^TAUO  >  **2  >  »*2 
ftETl.lRN 
CNO 

SUOROUTINC  RANC0F(XR.X|.AMP> 

THIS  function  returns  A  COMPLEX  RANfiOM  NUMBER.  (*R,Xl>.  WHERE  THE  RfPL 
AND  IMAOINARV  PARTS  ARE  INDEPENDENT  ZERO  MEAN  N0»'(MALLY  OlSTRieuTED 

random  variables  each  wiym  a  variance  of  .^*AMR**2. 

AROUMENTS 

XftsRETURNED  REAL  RANDOM  SAMPLE 
e  Xl*^RETfJRN£D  IM/VJtNARV  RANtiOrt  SAMPLE 

C  AMP’fSOL'AKg  root  CtF  TWO  TIMES  THE  VARt&NCE  OF  XR  AND  XI 
e  R£i?*.tlRED  functions 
C  60-RT  IFlOATINO  point  SCajAPf 

c  alch>  <  Float  tN»3  ptStNT  natural  LO<i> 

C  rand  (SAMPLE  C?F  UNIFORM  DISTRIBUTION  BETWEEN  JEHU  AND  0N£> 

C  COS  tPLOAriNO  POINT  COSINE! 

c  SIN  IFLOATINO  point  ■sINg! 

e 

AH-AMP*SOfi!T  t -AL oat  feAND*0»  »  » 

AN=ii  .  ’  a  3 1 *Ra.ND « 0  * 

XR=AM*cOStAN» 

Xt^AM*siNtANi 

RETURN 

END 

FUNCTICmM  RANDtl* 

e 

c  Tttis  ftWiCTICN  CifNERATES  PAfjDOM  NUMBERS  L*f{FC-i6!VV  DI-TfetDUTED 
e  lAPPwOplMATELY!  IN  twi  iNTERvtAL  Or  JERO  TO  ONE.  TwE  U^EH  SHOULD 
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C  RCPLACC  THIS  ALOOftlTHH  UITH  THE  BETTER  ONE(PRESUIAfiLT)  ON  HIS  OWN 
C  SVSTSN. 

C 

C  ARGUMENTS 
C  t-aiMMV  AROUMENT 
C  REOUIREO  COMMON  VARIABLES 

c  seeo-ranoom  number  SEEOi initial  value  must  be  000> 

C  RECJOIRED  FUNCTIONS 

C  OnOO  (DOUBLE  PRECISION  MOO  FUNCTION) 

C 

DOUBLE  PRECISION  SEED 
COMMiJN  /RAN/SEEO 

SEED«DMOO(4.0S:<3t2SD*O7«SEEO.2.6843S4S6O^> 

RAND- ( SEEO/ 2 . 6343S4570«Od ) 

END 

FUNCTION  P0UCR(0> 

C 

C  this  function  calculates  THE  GENERALIZED  POWER  SPECTRUM  AT  A  DELAY  CiF 
C  D  INTEORATEO  OVER  ALL  FREQUENCY. 

C 

C  ARGUMENTS 
C  0»DELAY{3EC» 

C  REO*.‘lR£D  common  VARIABLES 
C  CIH'ELAY  parameter 
c  FO*Fr<Ec»ucNcy  selective  bandwidth 

C  RECHJIREO  FUNCTIONS 
e  SORT  (FLOATING  POINT  SQUARE  ROOT) 

C  ABS  ( FLOAT  IN*:-  POINT  ABSiJLUTg  VALUE) 

c  EXP  (Floating  po)nt  natural  exponential) 
e 

COMMC»N  /PSD/ei.TA).tO.FO 
FPR lME=«FV*S'.'ftT  <  t .  *e  t  ♦•2  ) 

TPFPP'fc . 23S*FPR I ME *D 
?».70?I*(Ct-.TPFPD/ei) 

T«l .  / (  I .  * .  4?f;*4r»ABS <  J ) ) 

PT»t  ( .  74r5S^4*T'.0‘»SSTV-3l*T*.?480242)*T 
IFtf.LT.O.  XiO  TO  10 

PC'UER'^S.  I4)e*FMilMR*E»P4  =  .S*(TPPPD/e»  )**2)*PT 
RETURN 

to  P-^iUERsS.  »4|B*PPNlME*€XPt.5*Ct*»2-TPPPD)*{2.'FT*EtP(-2**i)i 

RETURN 
END 

FUNCTION  TC0N(T> 


Tl^tS  fUNCTIGN  calculates  TUg  COMPLEX  SIGNAL  TIHg  DECORRELATION 
PUNCTieN. 

ARGUMENTS 

T=TJM£  DI^PlA>-6m£NT(SEC» 

Rel^*J^Pfe  COMMON  L'AatABLES 

I a«Ar  DecGRRELATt^  TtPEiEEC'T 
nd^^lat  fade  pla>} 


SO 


C  NP-t.  F;.AT  Pf.JtNO 
C  NP.Or.t.  FftE'X'ESCV  SELECtlVE  FAOtN^ 

C  REOUJSED  FUNCTtONS 

C  E^f*  «FLCiATlNO  POINT  NflTURAl  EXPONENTJfit.1 
C  ASS  iFLOMtm  POINT  <^£OLU^E  VAUUE) 

c 

COWON  /PSO/C».TAyO.PO 

COWON  /OAT/NPP.OPP.NOP.KiP+NO.DS 

c 

C  CHECK  POP  PLAT  PAQinO 

c 

IP<NO.EQ. 1)00  TO  10 
TCOP-EXP(-<  f/TAU0>**2l 
return 

1C  TC0R»2. 146«A5S«T>/TA».I0 

TCOR«exP<-TCOft)*{l.*TCOA> 

RETURN 

END 

SUSROUTING  0£TD<TA|k..NT) 

C 

C  StfPROOTlNES  SETFn  PPT.  ANO  Pp&  fiCCO«PLliH  THE  PAST  FO'J^'lER  TRANSPORH 
C  ON  A  COMPLEX  ARRAV.C.  CR  AND  Cl  UTLL  t?€  USEO  A>  A  SHORT  HANC*  TO 
C  REPRESENT  THE  REAL  AND  IMAOINARV  PAATO  OP  C.  LET 
C  CR(I»  • 

c  cnt>  •  eu*T> 

C  DEFINE  THE  CORRESPONeeNCE  J^TUEEN  C.  CR.  AND  Cl.  THg  RESULT'S  ARE 
C  returned  in  TMfe  C  ARRAV.  TmE  R.JNSgR  OF  POINT'?  TO  DE  TRANSFORMED  f? 
e  NP»;t*NT.  THE  FA'iT  FO»>iiER  TRANSFORM  RE'.HJtRg-.  A  TAfcl.L  REFRiSENTgO  D¥ 

C  THE  COMPLEX  TADL  ARRAV,  MlNlKUH  LIMIT'S  FOR  THE  C  ANO  TAfiC  AS£ 

c 

C  DIMENSION  C<NP*NP>.TaftL(NP) 

c 

c  (CRisO.ciwo  at  the  Output  eouals  the  sum  oue»  i  from  t  to  np  of 
c 

C  teRU>*C'?S<2.*P{*M«M/NPl-Cl{ll*St««2.»Pl««*NMMPl.C6<  t>*S|N«i.*Pi*M.»N^ 
C  NP»*ei«l»*C00J2.*F'l*M*N.:'»jP»l 
€ 

C  UMCRE  THE  CCWPLE*  Ni_srdTltM.i  MAS  pEEN  USED.  tRgAL.  IKAOI.-^AA'iFI .  AND 
c  UHERE  N^-^l  UhEN  J.OE.  UAnD.  J.LT.f;P/?*l 
c  N^sO-NP-t  WHEN  ^*.0e  NP/2»l.AND.^.tf .NP 

c  UHEN  I.OE.I.ANFt. 

f.  H*:-NP-1  XSmEN  l.0£.NP/2*i.AN&t. 

e 

e  EUMi*!l».lT  INE  SETD  INI  ft  ALICES  THE  tA|i;,  AFlRAV  AND  NEED  PE  CAt4|&  OfiUV 
E  NHLN  NT  I?  C»AN>>ED.  's.»:^^OUt tN£i  EFT  AnD  P&O  UHgN  CAttED  SN 
C  ACCtWJH.tSH  THE  TfeA-NSFORM. 

e 

e  CALL  SEt&iTA&L.Nfl 

c 

f  CAt-L  PETiC.TAPL* 

e  €Atw  wiDtc* 

c 


St 


C  IN  THE  ABOVE  SUM.  M  IS  PROPORTIONAL  TO  THE  VARIABLE  REPRESENTED  BY  THE 
C  I  INDEX.  FOR  EXAMPLE.  IF  (CR< I ) .01 ( I ) )  IS  A  FREOUENCY  COEFFICIENT. 

C  THEN  ZERO  FREQUENCY  IS  REPRESENTED  BY  M  EQUALTO  ZERO  AND  I  EQUAL  TO 
C  ONE.  SIMILARLY.  N  IS  PROPORTIONAL  TO  THE  VARIABLE  REPRESENTED  BY  THE 
C  J  INDEX. 

C 

C  FAST  FOURIER  TRANSFORMS  ASSUME  CONTINUITY  OF  THE  FUNCTION  TO  BE 
C  TRANSFORMED  AND  THE  RESULT  AT  I (OR  J)  EQUAL  TO  NP/2+1  WHICH  REPRESENTS 

C  THE  EXTREMES  OF  M(OR  N).  THE  FUNCTIONS  ARE  ALSO  CONTINUOUS  BETWEEN 

C  (CRdl.CKD)  AND  (CR<NP)  .Cl  <NP) ). 

C 

C  REQUIRED  FUNCTIONS 

C  FLOAT  (FIXED  TO  FLOATING  POINT  CONVERSION) 

C  MOD  (FIXED  POINT  MOD) 

C  COS  (FLOATING  POINT  COSINE) 

C  SIN  (FLOATING  POINT  SINE) 

C 

C  THE  ARRAYS  IN  SETD.PFT.  AND  PDD  HAVE  BEEN  DIMENSIONED  TO  ONE 
C  INTERNALLY.  THIS  MAY  NOT  BE  ALLOWED  ON  THE  USER  SYSTEM.  IT  MAY  BE 
C  NECESSARY  TO  DIMENSION  TO  THE  LIMITS  IN  THE  MAIN  PROGRAM  OR  TO  USE 
C  VARIABLE  ARRAY  DECLARATORS. 

C  «««*»*»»•»*«•«*•***«*«*»'»•*•••*•*•*»««•«*«*»*««****««**«»•**«•***«•«** 

c 

DIMENSION  TABL(l) 

COMMON  /FFT/N2.NTI,NL 
DATA  PI/3, I4I592653S/ 

NTI«NT 

N2»2**NT 

NL=N2-1 

NN-N2/2-1 

TPI=>2.*PI/FL0AT(N2) 

TABL(1)=1. 

TABL(2)»0. 

DO  40  1=1. NN 
J»I  +  I 
II-O 

DO  20  N=1,NT 
II»II't-ll+M0D(J.2) 

20  J»J/2 

K-I+I+l 

TABL(K)=COS(TPI*FLOAT(lI ) ) 

40  TABL(K*l)=SlN(TPI*FLOAT(in) 

RETWN 

END 

CUBROUTINE  PFT(C.TABL) 

DIMENSION  Cin.TABLd) 

COMMON  /FFT/N2P,NTI.NL 

Nl=2 

N2oN2P 

DO  1  Iwl.NTl 

NN«-N2-N2 
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DO  2  J=2.N1,2 
NN=NN+N2+N2 
N3=NN+-N2 
DO  3  K=2,N2,2 
N4=»N3+K 

STR=C ( N4- 1 ) *TABL ( J-l) -C  <  N4 ) *TABL ( J ) 
STI=C(N4-1 )*TABU(J)+C<N4)*TABL(J-n 
N5=NN+K 

C(N4-1 )=C(N5-1 )-STR 
C(N4)=C(N5)-ST1 
C  <  N5- 1 >  =C ( N5- i ) +STR 
3  C<N5>=C<N5)+STI 

2  CONTINUE 
N1*N1+N1 

1  N2=N2/2 

RETURN 
END 

SUBROUTINE  PDD(C) 

C 

C  REQUIRED  FUNCTIONS 
C  MOD  (FIXED  POINT  MOD) 

C 

DIMENSION  C(l) 

COMMON  /FFT/N2.NTIfNL 
DO  1  1=1, NL 
J»I 
11=0 

DO  3  N=1.NTI 
II='II  +  II*M0D(J.2) 

3  J-U/2 
1F(II-I)1,1,7 

7  IN-I+I+I 

IIN-Il-*-II<.l 

0R=C<1N) 

DI-C(IN-*.l> 

C(IN)=C(IIN) 

C<IN+1)=C(IIN+1) 

C(nN)=DR 
C(HN+I  )=DI 
1  CONTINUE 

RETURN 
END 
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